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Abstract. We prove the so-called generalized Haff's law yielding the optimal alge- 
braic cooling rate of the temperature of a granular gas described by the homogeneous 
Boltzmann equation for inelastic interactions with non constant restitution coefficient. 
Our analysis is carried through a careful study of the infinite system of moments of the 
solution to the Boltzmann equation for granular gases and precise estimates in the self- 
similar variables. In the process, we generalize several results on the Boltzmann collision 
operator obtained recently for homogeneous granular gases with constant restitution co- 
efficient to a broader class of physical restitution coefficients that depend on the collision 
impact velocity. This generalization leads to the so-called -exponential tails theorem 
for this model. 



1. Introduction 

1.1. General setting. Rapid granular flows can be successfully described by the Boltz- 
mann equation conveniently modified to account for the energy dissipation due to the 
inelasticity of collisions. For such a description, one usually considers the collective dy- 
namics of inelastic hard-spheres interacting through binary collisions [ jTOl , p2| , The 



loss of mechanical energy due to collisions is characterized by the so-called normal resti- 
tution coefficient which quantifies the loss of relative normal velocity of a pair of colliding 
particles after the collision with respect to the impact velocity. Namely, if v and t>^ denote 
the velocities of two particles before they collide, their respective velocities v' and t>^ after 
collisions are such that 

{u ■ n) = —{u ■ n) e, (1.1) 

where the restitution coefficient e is such that ^ e ^ 1 and n G determines the 
impact direction, i.e. n stands for the unit vector that points from the f -particle center to 
the i;^-particle center at the instant of impact. Here above 

U = V — Vi,, u' = v' — f^, 

denote respectively the relative velocity before and after collision. The major part of the 
investigation, at the physical as well as the mathematical levels, has been devoted to the 
particular case of a constant normal restitution. However, as described in the monograph 
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[|rO|], it appears that a more relevant description of granular gases should deal with a 
variable restitution coefficient e(-) depending on the impact velocity, i.e. 

e := e(|M ■ n|). 

The most common model is the one corresponding to visco-elastic hard-spheres for which 
the restitution coefficient has been derived by Schwager and Poschel in [p2|]. For this 



peculiar model, e(-) admits the following representation as an infinite expansion series: 

oo 

e(|M.n|) = 1 + MeR^ neS^ (1.2) 

k=l 

where ^ for any k E N. We refer the reader to p2| ] for the physical considerations 
leading to the above expression (see also the Appendix A for several properties of e(-) in 
the case of visco-elastic hard-spheres). This is the principal example we have in mind for 
most of the results in the paper, though, as we shall see, our approach will cover more 
general cases including the one of constant restitution coefficient. 

In a kinetic framework, behavior of the granular flows is described, in the spatially situ- 
ation we shall consider here, by the so-called velocity distribution f(v, t) which represents 
the probability density of particles with velocity t; G at time t ^ 0. The time-evolution 
of the one-particle distribution function f{t,v),v E M'^, t > satisfies the following 

dtfit, v) = Qeif, m, v), fit = 0,v) = foiv) (1.3) 

where Qe{f, f) is the inelastic Boltzmann collision operator, expressing the effect of bi- 
nary collisions of particles. The collision operator Qe shares a common structure with the 
classical Boltzmann operator for elastic collision [|T1], ^ but is conveniently modified in 
order to take into account the inelastic character of the collision mechanism. In particular, 
Qe depends in a very strong and explicit way on the restitution coefficient e. Of course. 



for e = 1, one recovers the classical Boltzmann operator. We postpone to Section the 
precise expression of Qe- Due to the dissipation of kinetic energy during collisions, in the 
absence of external forces, the granular temperature 

£it)= [ fit,v)\v\'dv 



is continuously decreasing and is expected to go to zero as time goes to infinity, expressing 
the cooling of the granular gases. 

Determining the precise rate of decay to zero for the granular temperature is the main 
goal of the present work. The asymptotic behavior for the granular temperature was 
first explained in [|T5]] by P. K. Haff at the beginning of the 80's for the case of constant 
restitution coefficient, thus, it has become standard to refer to this behavior simply as 
Haff's law. 

The mathematical study of Boltzmann models for granular flows was first restricted to 
the so-called inelastic Maxwell molecules where the collision rate is independent of the 
relative velocity [|^, ^, H, [11], [1^]. Later, the mathematical investigation of hard- spheres 
interactions was initiated in [|14p for diffusively heated gases and continued in a series 
of papers []T7|, ITHI] where the first rigorous proof of the Haff's law was presented in the 



case of constant restitution coefficient. Additional relevant work in the existence and 
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Stability of the homogeneous cooling state can be found in [[T^, We refer to 
for a mathematical overview of the relevant questions addressed by the kinetic theory of 
granular gases and complete bibliography on the topic. 

From the mathematical viewpoint the literature on granular gases with variable resti- 
tution coefficient is rather limited. However, the Cauchy problem for the homogeneous 
inelastic Boltzmann equation has been studied in great detail and full generality in [fT7|], 
including the class of restitution coefficients that we are dealing with in this paper. For 
the inhomogeneous inelastic Boltzmann equation the literature is yet more scarce, in this 
respect we mention the work by one of the authors [jl]] that treats the Cauchy problem 
in the case of near-vacuum data. It is worthwhile mentioning that the scarcity of results 
regarding existence of solutions for the inhomogeneous case is explained by the lack of 
entropy estimates for the inelastic Boltzmann equation, thus, well known theories like the 
DiPerna-Lions renormaUzed solutions are no longer available. More complex behavior 
that involve boundaries, for instance clusters and Maxwell demons, a re well beyond of 
the present techniques. 

1.2. Main results and methodology. Physical considerations and careful dimensional 
analysis led P. K. Haff [ JT5| ] to predict that, for constant restitution coefficient, the temper- 
ature S{t) of a granular gas should cool down at a quadratic rate: 

S{t) = O (r^) as t ^ oo. 

Similar considerations led Schwager and Poschel [ p2| ] to conclude that, for the restitution 
coefficient associated to the visco-elastic hard-spheres (|1.2[), the decay should be slower 
than the one predicted by Haff, namely at an algebraic rate proportional to t^^^^. These 
considerations are precisely described in the main result of this paper where the key in- 
tuitive fact is that the decay rate of £(t) is completely determined by the behavior of the 
restitution coefficient e{\u ■ n\) for small impact velocity (Assumption (1) in |3.1|). Pre- 
cisely, our result is valid for restitution coefficient such that there exist some constants 
a > and 7^0 such that 

e{\u ■ n\) ^ 1 — a\u ■ n\'^ for |-u ■ n| ~ 

and reads as follows: 

Theorem 1.1. For any initial distribution velocity /o ^ satisfying the conditions given 
by ( |2.8| ) with /o G LP"(M^)/or some 1 < po < c>o, the solution f{t,v) to the associ- 
ated Boltzmann equation ([2.7|) satisfies the generalized Haff 's law for variable restitution 



coefficient e{-) fulfilling Assumptions and | 

c{l + ty^ ^£{t) ^C{l + t)-^, t^O (1.4) 

where £{t) = /j^a f(t, f )|f p df and c, C are positive constants depending only on e(-) 
and £{f)). 

We recover with Theorem [TTT] the optimal decay for constant restitution coefficient 



(7 = 0) given in nT9| ] and the one predicted for viscoelastic hard-spheres (7 = 1/5) in 



[E^ . The method of the proof has similarities to that of the constant restitution coefficient 



[]TP[] but technically more challenging. 
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The main tools to prove Theorem [Ll] are the following: 

The study of the moments of solutions to the Boltzmann equation using a generaliza- 
tion of the Povzner's lemma developed in [H]. 



• Precise estimates, in the same spirit of [|19|], of the solution to the Boltzmann equa- 
tion for p > 1. 

• For the previous item, the analysis is understood in an easiest way using rescaled solu- 
tions to ( [L3| ) of the form 

f{t,v) = V{tfg{T{t),V{t)v) 

where r(-) and V{-) are fixed time-scaling functions to be crafted depending upon 
the restitution coefficient. In the self-similar variables (r, w) the function gij^ w) is a 
solution of an evolution problem of the type 

9^5((r,w) + ^(r)V^ ■ iwg{T,w)) = Qe(r)i9,9) (1-5) 

for some ^(r) depending on the time scale r. The collision operator Qe{T){g, g) is as- 
sociated to a time-dependent restitution coefficient e(r) (see Section ^J]for details). In 
this respect we notice that one notable difference with respect to the case of a constant 
restitution coefficient treated in [ |T9| ] is that the rescaled collision operator depends on 
the (rescaled) time r, leading to a non-autonomous problem for g. This is the main 
reason why the construction of self-similar profile g independent of r obtained in [ [T^ ] 
(Homogeneous Cooling State) is not valid for non constant restitution coefficient. 

Let us explain in more details our method of proof: 

1 . We start proving in Sections 2 and 3 an upper bound for the decay of the energy. This 



shows that, for restitution coefficients satisfying |37I], the cooling of the temperature 
is at least algebraic. More precisely, under suitable assumptions on the restitution 
coefficient e(-), we exhibit a convex and increasing mapping ^'e such that 

^£{t) ^ -*e(^(t)) Vt ^ 0, 

at 

which leads to an upper bound for S{t) of the type 

£{t)^c{i + ty^ vt^o 

for some positive constant C > 0. 
2. The lower bound for the free cooling is much more intricate to establish and consists 
in proving that the cooling rate found above is optimal, i.e., there exists c > such 
that 

£{1)-^ c{l+t)~^i Vt^O. (1.6) 
A careful study of the moments of the solution to ([0|) shows that it suffices to prove 



a similar algebraic lower bound with some arbitrary rate, i.e. ( [Lq ) will hold if there 
exists A > and c > such that 

These two points are proved in the last part of Section 3. 
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3. To prove that the lower bound with some unprescribed rate A holds we use, as in 
[[r9|], precise estimates (p > 1) for solutions to (|T]3[) in self-similar variables. We 



craft a correct time scaling functions r(t) and V{t) such that ( [L^ ) is equivalent to 
0(r(t)) ^ c (here ©(•) denotes the second moment of g). Once this scale is fixed, 
the function g(T,w) satisfies the rescaled Boltzmann equation ( |1 .5^ with ^(r) — )■ 
as r — )■ oo. This is a major difference with the constant restitution coefficient case 
where ^(r) = 1. This technical difficulty is overcome proving that the L^-norms of 
5f(r) behaves at most polynomially with respect to r. For technical reasons which are 
peculiar to the inelastic interactions, noticed in [Q], we will restrict ourself to study 
L^-norms in the range p E [1,3). The details can b e found in Section 5. 

The derivation of precise estimates for the solution g(T, w) to ( |1.5| ) requires a careful 
study of the collision operator Qe and its regularity properties. We present in Section 
4 a full discussion of the regularity and integrability properties of the gain part of the 
collision operator ^ associated to a general collision kernel B(u, a) = ^{\u\)b{u ■ a) 
satisfying Grad's cut-off assumption (see Section 2 for definition). This Section is divided 
in five subsections starting with the Carleman representation of the gain operator Q^g. 
It is well-known ^ [T^] that such a representation is essential for the study of 



regularizing properties of the gain operator ^ when smooth assumptions are imposed 
on the kernel B{u, a). Our contribution in Sections 4.3 and 4.4 is to extend the existent 
theory to the inelastic case with variable restitution coefficient. Since the estimates of 
Section 4 will be applied for solutions written in self-similar variables, we make sure that 
such estimates are independent of the restitution coefficient. This allows us to overcome 
the technical problem of the time dependence of the gain operator in the self- similar 
variables. Additional convolution-like inequalities are derived in subsection 4.2 

assuming minimal regularity of the angular kernel 6(-). 

The final part of this work is devoted to the proof of propagation of exponential -tails 
where the full power of the Povzner's lemma is exploited. Much of the argument with a 
minor adaptation is taken from This important result is presented in the final Section 
for convenience and not because the machinery of Sections 4 and 5 is needed to prove it. 

Theorem 1.2 (L^ -exponential tails Theorem). Let B{u, a) = \u\b{u- a) be the collision 
kernel with b{-) satisfying ( |Z^ ) and b{-) G L'?(S^) for some q ^ 1. Assume that the 



variable restitution coefficient e(-) satisfies Assumptions |3. 1\ Furthermore, assume that 
/o satisfies (|2.8|), and that there exists tq > such that 



/ fo{v) exp (ro|t;|) dv < oo. 
Then, there exists some r ^ tq such that 

sup / f{t,v)exp{rV{t)\v\)dw<oo. (1.7) 



The function V(t) is the appropriate scaling, depending solely on the restitution coeffi- 
cient, given in ( p.l6D . 
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1.3. Notations. Let us introduce the notations we shall use in the sequel. Throughout the 
paper we shall use the notation (•) = a/1 + | ■ p. We denote, for any G M, the Banach 



space 



Ll=\f-- 



measurable : 



\f{v)\ {v)''dv <+oo\. 



More generally we define the weighted Lebesgue space L^(]R^) {p E [1, +oo), i] eW) 
by the norm 



\fiv)f{vrdv 



i/p 



1 ^ p < oo 



while ||/||lc»(k3) = ess - sup^^^3\f{v)\{v)'i forp = oo. 

For any A; G N, we denote by H'' = H''{R^) the usual Sobolev space defined by the 
norm 

i/p 

where dl denotes the partial derivative associated with the multi-index j E N^. Moreover 
this definition can be extended to H'^ for any s ^ by using the Fourier transform J^. The 
binomial coefficients for non-integer p ^ and /c G N are defined as 



p 

vJ iil2 



p{p - 1) . . .{p - k + 1) 



A; > 1, 



1. 



2. Preliminaries 

2.1. The kinetic model. We assume the granular particles to be perfectly smooth hard- 
spheres of mass m = 1 performing inelastic collisions. Recall that, as explained in the 
Introduction, the inelasticity of the collision mechanism is characterized by a single pa- 
rameter, namely the coefficient of normal restitution ^ e ^ 1 which we assume to be 
non constant. More precisely, let {v,Vi,) denote the velocities of two particles before they 
collide. Their respective velocities after collisions v' and are given, in virtue of ( |1.1| ) 
and the conservation of momentum, by 

1 + e ^ ^ 1 + e ^ ^ 
v' = V — {u ■ n)n, v'_^ = v^, -\ — {u ■ n)n, (2.1) 

where the symbol u stands for the relative velocity u = v — v^, and n is the impact 
direction. From the physical viewpoint, a common approximation consists in choosing e 
as a suitable function of the impact velocity, i.e. e := e{\u ■ n\). The main assumptions 
on the function e(-) are listed in the following (see \^y. 

Assumptions 2.1. Assume the following hold: 

(1) The mapping r E M+ e(r) G (0, 1] is absolutely continuous. 

(2) The mapping r E M+ — )■ 'i9(r) := r e(r) is strictly increasing. 
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Figure 1 . Restitution coefficient for viscoelastic hard-spheres given by 
Eq. Q witha = 0.12. 



Further assumptions on the function e(-) shall be needed later on. Given assumption 
(2), the Jacobian of the transformation ( jOI ) can be computed as 



J :-- 



d{v, 



^ ^ de ^ di9 ^ 

\u -nl + \u ■ n\ — (\u ■n\) = — (Im ■ n|) > 0. 
dr dr 



In practical situations, the restitution coefficient e(-) is usually chosen among the follow- 
ing three examples: 

Example 2.2 (Constant restitution coefficient). The most documented example in the 
literature is the one in which 

e(r) = Co G (0, 1] for any r ^ 0. 

Example 2.3 (Monotone decreasing). A second example of interest is the one in which 
the restitution coefficient e(-) is a monotone decreasing function: 

e(r) = Vr ^ (2.2) 

where a > 0, > are two given constants. 



Example 2.4 (Viscoelastic hard-spheres). This is the most physically relevant model 
treated in this work. For such a model, the properties of the restitution coefficient have 
been derived in [[T^, ^ where representation ( |1.2D is given. It also accepts the implicit 
representation 



(2.3) 



where a > Ois a suitable positive constant depending on the material viscosity ( see Figure 
1). 
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In the sequel, it shall be more convenient to use the following equivalent parametriza- 
tion of the post-coUisional velocities. For distinct velocities v and v^,, let m = |^ be the 
relative velocity unit vector. The change of variables 

a = u — 2{u ■ n)n E 

provides an alternative parametrization of the unit sphere §^ for which the impact velocity 
reads 



1 — u ■ a 

\u ■ n\ = \u\ \u ■ n\ = \u\\l . 

2 

Then, the post-coUisional velocities {v', v'^) given in (pTTj) are transformed to 



, — \u\a , — \u\a ,^ ,^ 

v' = v- {3 v[ = v^ + P ^ (2.4) 



where 



/3 = /3(HV^)=^e(ii]- 



In this representation, the weak formulation of the Boltzmann collision operator QB,e 
given a collision kernel B(u, a) reads 

QBAf^9){v)ij{v)dv =]- I f{v)g{v^)AB,e['^]{'^^v^) (2.5) 

for any suitable test function tp = ^p{v). Here 



with v', v'^ defined in (|Z4|). We assume that the collision kernel B{u, o) takes the form 

B[u,a) = <^{\u\)b{u ■ a) 

where $ ( • ) is a suitable nonnegative function known as potential, while the angular kernel 
b{-) is usually assumed belonging to (— 1, 1). For any fixed vector u, the angular kernel 
defines a measure on the sphere through the mapping cr G i— )■ b(u ■ a) E [0, oo] that we 
assume to satisfy the renormalized Grad's cut-ojf hypothesis 

ll&llLi(s^) = 27r||6|L.(.i,i) = l. (2.6) 

The most relevant model in our case is hard-spheres which correspond to $(|m|) = \u\ and 
h{u- a) = We shall consider the generalized hard-spheres collision kernel for which 
$(|m|) = |m| and the angular kernel is non necessarily constant but satisfying ([2.6D. For 
the particular model of hard-spheres interactions, we simply denote the collision operator 

QB,e by Qe. 
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2.2. On the Cauchy problem. We consider the following homogeneous Boltzmann equa 
tion 

fio,v) =foiv), vem:' 

where the initial datum /o is a nonnegative velocity function such that 

Uv)dv = l, [ fo{v)vdv = and / fo{v)\v\'' dv < oo. (2.8) 



There is no loss of generality in assuming the two first moments conditions in d^q ) due to 
scaling and translational arguments. We say that a nonnegative f = f{t, v) is a solution 

to if / G C([0, oo), Ll(R^)) and 

dt [ (f{t,v)dMt,v) + QBMJ)it,^)Ht,v)]dv=[ fo{v)ij{0,v)dv 

holds for any compactly supported G ( [0, oo) x M'^) . Under the Assumptions PTT] , the 
assumptions HI and H2 of [[T^] are fulfilled (with the terminology of [[1^], we are deal- 
ing with a non-coupled collision rate and, more precisely, with the so-called generalized 
visco-elastic model, see [jlEp, p. 661). In particular, Theorem 1.2] applies direclty 
and allows us to state: 



Theorem 2.5 (Mischler et a/.). For any nonnegative velocity function /o satisfying ( |2.8| ), 
there is a unique solution f = f{t, v) to ( p.?] ). Moreover, 

I f{t,v) dv = l, I f{t,v)vdv = Vt^O. (2.9) 

2.3. Self-similar variables. Let us discuss precisely the rescaling using self-similar vari- 
ables. Let f{t, v) be the solution to ( pJ] ) associated to some initial datum /q satisfying 



( |2.8| ) and collision kernel 

B{u,a) = ^{\u\)b{u-a) 
with 6(-) satisfying (^^). The rescaled solution g = g{T, w) is defined such that 

f{t,v) = V{tfg{T{t),V{t)v) (2.10) 

where r(-) and V{-) are time-scaling functions to be determined solely on the behavior 
of the restitution coefficient in the low impact velocity region. Since these are scaling 
functions they are increasing and satisfy r(0) = and ^(0) = 1. One has 

1= / f(t,v)dv= / g{T{t),w)dw Vt ^ 

and g{0, w) = fo{w). Furthermore, some elementary calculations show that the function 
g{T, w) satisfies 



V{t)-'Qe{f,f){t,v) = f{t)V{t)drg{r,w) + ■ {wg{T,w)) 



(2.11) 



w=V{t)v 
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where the dot symbol denotes the derivative with respect to t. Moreover, the expression 
of the collision operator in the self-similar variables is 



where the rescaled collision kernel Br is given by 



The rescaled restitution coefficient has been defined by 

Cr : (r,t) ^ e,(i)(r) := e [y^ for r ^ 0, 0. 

Since the mapping t e IR+ i — > T(t) G is injective with inverse C, one can rewrite 
equation ( |2.1 ID in terms of r only. Thus, g{T, w) is a solution to the following rescaled 
Boltzmann equation: 

\{r)drg{T,w) + i{T)V^ ■ {wg{T,w)) = QBrM9,9)ir,w) Vr > (2.12) 

with 

A(-) = r(C(-))nC(-)) and ^O) = ^CH), 
and model parameters 

BAu,a) = Viar))^{yl^)Ku.a) and ^rir) = e ^^^y (2.13) 

Notice that, for generalized hard-spheres interactions (i.e. whenever = one 

has Br = B. For true hard-spheres interactions, i.e. = one simply denotes the 
rescaled collision operator by Qe^ • In addition, observe that the rescaled operator depends 
on time, and therefore, g is a solution to a non-autonomous problem. 

2.4. Povzner-type inequalities. We extend in this section the results of [Q] and [[T^p to 
the case of variable restitution coefficient satisfying pTT] . We consider a collision kernel of 
the form 

B{u,a) = <^{\u\)b{u-a), 



with angular kernel b{-) satisfying the renormalized Grad's cut-off assumption (|2.6D. Let 
/ be a nonnegative function satisfying ( |2.9| ) and il){v) = p) be a given test-function 



Qb A f J) iv)Hv)dv = I I f{v)f{v,)AB,e[i^]iv,v,) dv. dv 

R^xn 



with \1' nondecreasing and convex. Then, Eq. ( p31 ) leads to 

1 
2 

AsAi^Kv^v^) = HH) {A+BA'^]iv,v.) - A-^A'^]iv,v.)) 
AtAmv,v.) = [ {n\vf) + n\v'f))b{u-a)da. 



with 
where 
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Using ( [2 .61 ) we also have, 



Following [pp, we define the velocity of the center of mass U 



so that 



\u\ \u\ 
v' = U + ^u, vl = U - ^uj with w = (1 - /3)u + /3a. 

Recall that for any vector x E M^, we set x = -pr . When e, or equivalently f3, is constant, 

\x\ 

the strategy of consists, roughly speaking, in performing a suitable change of unknown 
cr — )■ cD to carefully estimate e[^]- variable /3, such strategy does not apply directly. 
Instead, observe that ^ 1 and, since ^ is increasing, one has 

^(|t;f) + ^(|<n ^ ^ ( \u\^ + + \u\\U\U ■ tu] +^ (\U\^+^-^-\u\\U\U-u 



4 



* I E I + ^ ^ 



4 



where we have set E := \v\^ + \v^\^ = 2|f/p + and ^ = 2^^^. Since ^(■) is convex 
the mapping 

■q)Q(t) = ^(x + ty) + ^(x - ty) 

is even and nondecreasing for t ^ and x,y E M (see Therefore, using that ^ ^ 1 
one gets 



(2.14) 



In the case that U ■ a > it follows that 



U-uj 



[1- /3)U -11 + (3U -a ^{1- (3) + i3U -a. 



thus, using the fact that \l'o(i) is even and nondecreasing for t ^ 0, we conclude from 
(Og)that 



^{\v'\^) + ^{\vX) [ E — 1+^1^- 



When U ■ a ^0 a similar argument shows that 



^{\v'\^) + '^{\v'^\^) \ E — \+^\E- ^ 
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Hence, setting b(s) = b{s) + b{—s) and using these last two estimates with the change of 
variables a — )■ — a we get 























y{C/-cr>0} 





^ E 



(3-I3U -a 



h{u ■ a) da 



,3 + t/- a 



1-U-a 



h{u ■ a) dcr, 



(2.15) 



where the second inequality can be shown writing 



and, 



(3-l3U-a 1 /I 13 



l-U 



a 



The term in parenthesis is maximized when (3 = 1/2, thus the monotonicity of \l'o implies 
the result. 

Next, we particularize the previous estimates to the important case '^{x) = x^. This 
choice will lead to the study of the moments of solutions. 

Lemma 2.6. Let q ^ 1 be such that b G L''(S^). Then, for any restitution coefficient e(-) 
satisfying Assumptions |2.i| and any real p ^ 1, there exists an explicit constant Hp > 
such that 



]. (2.16) 



+ ^pm' + \v.\'Y-\v\'^-\v4'r 

This constant Hp has the following properties: 

(1) «:i ^ 1. 

(2) For p ^ 1 the map p ^ Hp is strictly decreasing. In particular, Kp < Iforp > 1. 

(3) Kp = O [l/p^/'^') for large p, where l/q + = 1. 

(4) For q = 1, one still has Kp \ as p ^ oo. 

Proof. Let '^p{x) = x'^. From ( |2.15D , one sees that 

AlA^p]iv,v^)^KpE^ 
where we recall that = If P + If^P and we set 



Kp = sup 



3 + U -a 

4 



1 - U 



a 



4 



b(u ■ cr) dcr. 



(2.17) 



U,u JU-a^O 

It is clear that the above inequality yields ( |2.16D . Let us prove that Kp satisfies the afore- 
mentioned conditions. First, we use Holder inequality to obtain 



Kp ^ 4:71 



1 r 



-1 



3 + s 



1 - s 



l/q' 



ds 



IGtt 



< 



Il'j(S2) 



{q'p+iy/i' 
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This proves that Kp is finite and also yields item (3) for g > 1. For items (1) and (2) 
observe that the integral in the right-hand-side ( [2.15D is continuous in the vectors U, u G 
This can be shown by changing the integral to polar coordinates. Thus, the supremum 
in these arguments is achieved. Therefore, there exist Uq, uq E §^ (depending on the 
angular kernel b) such that 



/- 



3 + Up - a 
4 



1-Uo-cj 



b{uo ■ a) da. 



A simple computation with this estimate shows that ki 



IL1(§2) 



1. Moreover, the 



integrand is a.e. strictly decreasing as p increases and this proves (2). Finally, let p — )• oo 
in this expression and use Dominated convergence to conclude (4) for the case q = 1. □ 

The above lemma is the analogous of [|^, Corollary 1] for variable restitution coefficient 
e(-) and it proves that the subsequent results of extend readily to variable restitution 
coefficient. In particular, Lemma 3] readsQ: 

Proposition 2.7. Let f be a nonnegative function satisfying (|2.9|). For any p ^ 1, we set 



rrir, 



f{v)\v\'^fdv. 



Assume that the collision kernel B{u, a) = \u\b{u ■ a) is such that b{-) satisfies ( [2.6[ ) with 
b{-) e L'^{S>^) for some q ^ I. For any restitution coefficient e(-) satisfying Assumptions 
|2.i| and any real p ^ 1, one has 

[ QBAfJ')iv)M'''dv^-{l-Kp)mp+y2 + ^pSp, (2.18) 



where, 



^^^] denoting the integer part of 2±i and Kp being the constant of Lemma 



E 

k=l 



Inequality ( |2.18D was introduced in because the term Sp involves only moments of 
order p — 1/2. Thus, the above estimate has important consequences on the propagation 
of moments for the solution to (p7|) (see Section p| for more discussion). 



3. Free cooling of granular gases: generalized Haff's law 

We investigate in this section the so-called generalized Haff's law for granular gases 
with variable restitution coefficient. More precisely, we aim to derive the exact rate of 
decay of the temperature £{t) of the solution to Eq. In this section, we exclusively 

study the generalized hard-spheres collision kernel. 

B{u, cr) = \u\b(u ■ a) 



'Notice that, though stated for hard-spheres interactions only, Lemma 3] applies to our situation 
thanks to the above Lemma 2.6 and Lemma 1]. 
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where b{-) satisfies ( |2.6| ) but generalization to the so-called variable hard-spheres interac- 
tions (i.e. = \u\'^ for s ^ 0) is easy to handle. Let /o be a nonnegative velocity 
distribution satisfying ( ^l8| ) and let f{t, v) be the associated solution to the Cauchy prob- 



lem (|2.7[). We denote its temperature by S{t), 

£it)= [ fit,v)\v\'dv. 

The conditions ( |2.8| ) implies that sup^^jQ < oo. Indeed, the evolution of £{t) is gov- 
erned by 



\u\ X 



— |t>*|^ J b{u ■ a) da dVi, dv 



i\2 I I / |2 



where we applied ( |2.5| ) with ) = |t>| . One checks readily that 



1 1 — M ■ 0" 



— \u\ 



1 -eM kl 



1 — u ■ a 



so that 



/" 1 — M ■ cr 




4 


;■-( 







1 — M • (7 



6(u ■ a) da. 



We compute this last integral over §^ (for fixed v and i;^,) using polar coordinates to get 



\u\ 



f I — u ■ a 









1 — u ■ a 



h{u ■ a) da 



27r|M|3 / (1 - eWu\y)) 6(1 - 2y^)y^ dy = 

^0 



u\ 



where we have defined 



Vr > 0. 



(3.1) 



*e(r) := 2Txr^'^ [ (1 - e{^/^zf) b {l - 2z^) dz, 
Jo 

la other words, the evolution of the temperature £(t) is given by 

4-£it) = - I f{t,v)f{t,v,)^,{\u\^)dvdv, < 0, t ^ 0. 

JlR3xR3 

Li addition to Assumptions pTTI , we assume in the rest of the paper that the restitution 
coefficient e(-) satisfies the following: 

Assumptions 3.1. Assume that the mapping r e(r) G (0, 1] satisfies Assumptions |Z7| 
and 
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(1) there exist a > and 7^0 such that 

e(r) ^ 1 — ar'^ for r ~ 0, 

(2) liminfr^ooe(r) = cq < 1, 

(3) b{-) e Li{S'^) for some q ^ 1, and 

(4) the function r > i — y ^'e('") defined in ( |3.1[ ) is strictly increasing and convex 
over (0, +oo). 

Remark 3.2. For hard-spheres interactions, b(u ■ cr) = j-, thus, reduces to 



{l-e{y f) y'dy, 



r > 0. 



We prove in the Appendix that Assumptions ^J] are satisfied for the viscoelastic hard- 
spheres of Example with 7 = 1/5. More generally, in the case of hard-spheres inter- 
actions, assumption (4) is fulfilled if e{-) is continuously decreasing (see Lemma A. in 
Appendix A). For constant restitution coefficient e{r) = cq, these assumptions are trivially 
satisfied. 

3.1. Upper bound for £{t). We first prove the first half of Haff's law, namely, the tem- 
perature £{t) has at least algebraic decay. 

Proposition 3.3. Let fo be a nonnegative velocity distribution satisfying (|2.8|) and let 



f{t, v) be the associated solution to the Cauchy problem ( |2.7D where the variable restitu- 
tion coefficient satisfies Assumptions pJ} Then, 

at 

Moreover, there exists C > such that 

S{t) ^C{l + t)~^ Vt^O. 
Proof Recall that the evolution of the temperature is given by 

—£{t) = - / f{t,v)f{t,v,)^,{\u\^)dvdv., t ^ 0, 



(3.2) 



(3.3) 



where u = v — v^. Since ^e(| ■ p) is convex according to Assumption |3TT| (2) and 
/(t, Vi,) dVi, is a probability measure over MP, Jensen's inequality implies 

2^ 



fit,v,)^M^)dv,^^, 



V 



v^f{t, v^) dv^ 



where we used ([Z^. Applying Jensen's inequality again we obtain 
and therefore. 



dt 



£{t) ^ -*e(^(t)) Vt ^ 0. 
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Note that ^'e(") is strictly increasing with limj._>o ^ei^) = 0, this ensures that 

hm S{t) = 0. 

t—^co 

Moreover, according to Assumptions pTI] (i j, it is clear from (p.ip that 



^e(x) ~ CyX 2^ for X ~ 0, 



where the constant can be taken as = 2-na y^+^b{l-2y^) dy < oo. Since ^(t) ^ 0, 



there exists to > such that ^'e(£^(t)) ^ iC^£^(t)^2^ for all t ^ to which implies that 

dt ^ ' 2 ^ ' " 



This proves ( p^ D 



Example 3.4. In the case of constant restitution coefficient e(r) = cq E (0, 1) for any 
r ^ 0, /or hard-spheres interactions, one has 

o 

TTiMi', one recovers from ( |3.2| ) decay of the temperature established from physical 
considerations (dimension analysis) in [ JT5| ] and proved in [JTP|], namely, £{t) ^ C(l+t)~^ 
/or Zar^e t. 

Example 3.5. For f/ie restitution coefficient e(-) associated to viscoelastic hard-spheres 



(see Example ^.4\), one has 7 = 1/5, ^/ze above estimate ( p.2| ) /eaJ* to a Jecay o/?/ze 
temperature faster than (1 + t)^^/^ which is the one obtained in [ p2| | f^ee a/50 [|T0|] j/rom 
physical considerations and dimensional analysis. 



Notice that, since £{t) —i- as t —t- 00, it is possible to resume the arguments of [ |T8| , 
Prop. 5.1] to prove that the solution /(t, v) to ( |2.7D converges to a Dirac mass as t goes to 
infinity, namely 



(j„_n wcais-iy ^ 111 ivi [ii^'^'^ 

t—^ca 

where A'/^(R^) denotes the space of normalized probability measures on M^. We shall not 
investigate further on the question of long time asymptotic behavior of the distribution 
/(t, v) but rather try to capture the very precise rate of convergence of the temperature to 
zero. 



/(t, v) — y 6^=0 weakly * in M\M.^ 



Using the Povzner-like estimate of Section it is possible, from the decay in £{t), to 
deduce the decay of any moments of /. Indeed, for any t ^ and any p ^ 1 we define 
the p-moment of / as 



mp{t) := / f{t,v)\v\^Pdv. (3.4) 
Corollary 3.6. Let /o be a nonnegative velocity distribution satisfying ( |2.8| ) and let f{t, v) 



be the associated solution to the Cauchy problem ( |2.7| ) where the variable restitution 
coefficient satisfies Assumptions For any p^ 1, there exists Kp > such that 

mp{t) ^ Kp{l + t)~^ Vt^O. (3.5) 
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_ 2 

Proof. Set u{t) = (1 +t) . We prove that, for any p ^ 1, there exists Kp > such that 
mp{t) ^ KpU^if) for any t ^ 0. Observe that using classical interpolation, it suffices to 
prove this for any p such that 2p G N. We argue by induction. It is clear from Proposition 
?3| that estimate (p3|) holds for p = \. Let p > 1, with 2p E N, be fixed and assume that 



for any integer 1 ^ j ^ p — 1/2 there exists Kj > such that mj{t) ^ KjU^{t) holds. 
According to Proposition ^77| 

-nmp{t)= / QB,e(/,/)(t,t^)|t^|'^dt; ^ -(l-Kp)mp+i/2(t) + «p5p(t), (3.6) 

where 



r p+i 
2 



'S'pW = X] ( A; ) ("^'=+i/2(^) ^p-k{t) +mk{t) mp_fc+i/2(t)) , Vt ^ 0. 
fc=i ^ ^ 

For p ^ 2, the above expression S'p(i) involves moments of order less than p — 1/2. The 
case p = 3/2 is treated independently. 

• Step 1 (p = 3/2). In this case ( |X6[ ) reads 

^m3/2(t) ^ -(1 - f^3/2)m2it) + m3/2(t)mi/2(t) + ^'(t) Vt ^ 0. (3.7) 
Let Khe a positive number to be chosen later and define 

f/3/2(t) ■.= my2{t)-Ku{tf^. 

Using ( p.7D one has 



^ _(1 _ Ky2)m2{t) + m3/2(t)mi/2(t) + S\t) + -^(1 + ty- 
at 1 + 7 

From Holder's inequality. 



4+7 
+ 7 . 



4+7 
1+7 



m3/2(t) ^ ^/S{t)Vm2{t) and mi/2(t) ^ ^/£it) Vt ^ (3.8) 

hence, 

df/3/2 '"A/ii't) > o 

-^{t) ^ -(1 - '^3/2)^|y^ + v^m3/2(t) + S\t) + + t)-^ 

_ 2 

Since < C(l + t) 1+^ , there exist a, 6, c> such that 

^^(t) ^ -am^/2(t)(l + t)TTT + 6 (1 + t)-TTT 

+ c(l + t)"^m3/2(t) + -^i^(l+t)"^ Vt>0. (3.9) 

1 + 7 

Inequality ( [J!9| ) implies the result for the case p = 3/2 provided is large enough. 
Indeed, choose K so that m3/2(0) < Ku^/'^{Q) = K. Then, by time-continuity of the 
moments, the result follows at least for some finite time. Assume that there exists a time 
> such that m3/2(t^) = Ku^/^{Q = K{l + Q~^, then implies 



{Q ^ ( -aK' + 6 + cir + K ] (1 + < 
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whenever K is large enough. Thus, ( |33| ) holds for p = 3/2 choosing '■= K. 

• Step 2 (p ^ 2). The induction hypothesis implies that there exists a constant Cp > 
such that 

where Cp can be taken as 

Cp = ^ i ^ j [Kk+1/2 Kp^k + Kk Kp^k+1/2) ■ 
k=i ^ ^ 

Furthermore, according to Jensen's inequality mp+i/2(t) ^ mp^^^^^(t), for any t ^ 0. 
Thus, from (|3.6|), we conclude that 

-^rupit) ^ -(1 - Kp)mJ+i/2p(^) ^ CpM(t)P+^/2 Vt ^ 0. 

(J.6 

Arguing as in Step 1, for some i^' > to be chosen later, we define 

Up{t) := mp{t) - Kuitf. 

In this way, 

dr i + 7 

Then, if is such that f/p(0) < 0, the result holds at least for some finite time. For any 
> such that Up{t^) =0, one notices then that 

1+7 < 



Up{U) ^ (^-(1 - Kp)K^^^v +KpCp + (1 + U 



provided K is large enough. This proves (^3J) for any p ^ 1. □ 

3.2. Lower bound for preliminary considerations. The next goal is to complete 
the proof of Haff's law by showing that the cooling rate ( |3l2[ ) is optimal under Assump- 
tions IsTTI. Thus, we have to show that there exists C > such that 



£{t)^c{i + ty— vt^o. 

First, we prove the following result that simplifies our endeavor. 

Theorem 3.7. Assume a non constant (7 > Oj restitution coefficient e(-) satisfying As- 
sumptions If there exist Co > and A > such that 

S{t)^Co{l+t)-^ Vt^O, (3.10) 

then there exists Cp > such that 

mp{t) ^ CpSP{t) for any t ^ and p ^ 1. (3.11) 

As a consequence, there exists C > such that 

£{t)^C{l+ty^ Vt^O. (3.12) 
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Proof. According to Assumption |37I] (1) 

3+7 

^e{x) — C^x~ for X ~ 0. 
In addition, Assumption ^7T] (2 j implies that there exists C5 > such that 

^e(a;) ^ ChX^^"^ for large X, 

where the constant can be taken as Ch = 27r(l — Cq) 6(1 — 2z'^)z^ dz. Thus, there exists 
another constant C > such that 

*e(2;) ^ Cx— Wx > 0. (3.13) 

Then, from ( P3| ) one deduces that for any e > and p > 



< C (^£^m3/2(t) + + ^)"^) Vt ^ 0. 



where we have used Corollary |3]^ for the second inequality. In particular, using ( |3.1(]| ) 
and the fact that £{t) is a non increasing function, one can choose p sufficiently large so 
that 

for some positive constant Cp. In other words, for any 5 > there exists > such that 

- ^ '5^3/2(t) + Cs£{tf/^ Vt ^ 0. (3.14) 

at 

With this preliminary observation, the proof of ( p.l ID is a direct adaptation of that of 
Corollary p^ . Here again, by simple interpolation, it is enough to prove the result for any 
p such that 2p G N and argue using induction. The result is clearly true for p = 1 with 
Ci = 1. For p = 3/2, let /"T > be a constant chosen later and define 

u^/2{t)=m^/2{t)-KS{tf\ 

Thus, from (|3]3) 

d 3 ci 

^^^3/2(t) ^ -(1 - ^^■i/2)m2{t) + m3/2(t)mi/2(t) + S\t) - -K^/S{i)—S{t). 

Using ( p.8p one deduces from ( [3.14p that, for any 6 > there exists > such that 

^«3/2(t) ^ -(l-'^3/2)^^||y^+(l + Ik6^ m3/2(t)V^+(l + IkCs^ £\t). 

Fix 5 = ^-^^ and choose K > such that ^3/2(0) < 0. If > is such that M3/2(t*) = 
0, then the following holds 

^«3/2(t.) ^ (-i^i^^ +i^K + l + IkCs^ ) ^(t.)^ < 0, 
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provided K is sufficiently large. This proves ( |3.1 ID for p = 3/2 with C3/2 := K. The 
case p ^ 2 follows in the same lines of the proof of Corollary ^]6| interchanging the roles 
of S{t) andn(t). 

To conclude the proof, observe that according to ( |3.13| ) and (P3|), there exists C > 
such that 

Then, applying ( |3.11| ) with p = one deduces that there is > such that 

A simple integration of this inequality yields ( [3.12[ ). 

Remark 3.8. For constant restitution coefficient, e 
hold anymore. However, for some Ce > Owe have 



Vt ^ 0. 



□ 



Co, since 7 = 0, (3.14) does not 



~—£it) ^ C,m3/2{t) Vt ^ 0. 

Assuming that Cq — 1 (quasi-elastic regime) the constant Ce is small, thus, the argument 
above can be reproduced to prove that the conclusion ofProposition^^still holds. Recall 
that for •y = the second part ofHaff's law ( |3.12| ) has been proved in [|T9|, Theorem 1 .2]. 

In order to prove that ( |3.10| ) is satisfied for some A > 0, we will need precise 
estimates, following the spirit of [JT9|], for the rescaled function g given in Section [23 



The idea to craft the correct time- scaling functions r(-) and V(-) is to choose them such 
that the corresponding temperature of g is bounded away from zero. Indeed, for any 
r > 0, define 

g{T,w)\w\'^ dw. 



Since, 
we choose 



@(t) 



£{t) = V{t)-^&{T{t)) vt ^ 0, 



V{t) = (1 + t) — 



Vt > 0. 



(3.15) 



(3.16) 



In this way, ( |3.12| ) is equivalent to ©(r(t)) ^ C for any t ^ 0. Notice that ( pl2| ) immedi- 
ately translates into 

sup e(r(t)) < 00. (3.17) 



Moreover, for simplicity we pick r(t) such that r(t)V{t) = 1, therefore for 7 > 0, 

"* ds 7 + 1 



r(t) 



7 

1+7 



(3.18) 



V{s) 7 

which is an acceptable time-scaling function. Thus, the rescaled solution g{T, w) satisfies 
( pl| ) with A(r) = 1, 



7r + (1 + 7) 



and er{r) = e r ( 1 + 



7 



7 + 1 



-r 



-1/7^ 



(3.19) 
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If 7 = the restitution coefficient is constant [fT9[], in particular = e, and the rescale 



reads V{t) = 1 + t and r(t) = hi(l + t). In such a case, ^(r) = 1. 

To complete the proof of Haff 's law, one has to perform a careful study of the properties 
of the collision operator Qe in Sobolev or spaces 1 < p ^ oo. 

4. Regularity properties oe the collision operator 



In this section the regularity properties studied originally for the elastic case in Qlq , 



2T| , p5| ] and later for the constant restitution coefficient in [ |T9| ] are generalized to cover 
variable restitution coefficients depending on the impact velocity. The path that we follow 
closely follows 

4.1. Carleman representation. We establish here a technical representation of the gain 
term ^ which is reminiscent of the classical Carleman representation in the elastic 
case. More precisely, let B{u, a) be a collision kernel of the form 

B{u,a) = <^{\u\)biu ■ a) 

where $(■) ^ and b(-) ^ satisfies (|2^). For any = ip{v), define the following linear 
operators 

S±{tlj){u)= / tlj{u^)b{u ■ a) da, Vu G M^ (4.1) 



the symbolos u and m+ are defined by 



u := p \ \u \\l 1 , and u .= u — u . 

Lemma 4.1. For any continuous functions ip and (f, 

ip{u)S-{'ilj){u)^{\u\) du = / ilj{x)TBi^)ix)dx 
where the linear operator T b is given by 

Tb{<^){x)= / B{z + a{r)u,a{r))ip{a{r)uj + z)dTXz, 

X = ru, r ^0, u e (4.2) 

Here dir^ is the Lebesgue measure in the hyperplane uj-^ perpendicular to u and a{-) is 
the inverse of the mapping s H- sf3{s). Moreover, 



with ■&{■) defined in As sumption ^J\( 2) and ■&'{■) denoting its derivative. 
Proof For simplicity assume that $ = 1. Define 

I .= ip{u)S^{%jj){u) du = / (p{u) du / ip{u~)b{u ■ a) da. 



§2 
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For fixed m e M^, we perfomi the integration over using the formula 

F{ da = — / 5{\x\^ - x-u)F{x)dx 

2 V 2 y \u\ 

valid for any given function F. Then, 

V9(-u)|'u|~"'"(5(|xp — X ■ u)ip(^xf3{\x\))b fl — 2-j— L j dxdu. 



Setting now m = z + x we get 

I = A j Lp{x + z)\x + z\~^5{x ■ z)il)[x(3{\x\))h (l - 2 \ dzdx. 

Jr3xR3 V \X + z\ J 

Keeping x fixed, we remove the Dirac mass using to the identity 

f F{z)5{x ■ z)dz = [ F{z)d7i,, 
Jr3 \x\ Jx± 

which leads to 

Perform the x-integral using polar coordinates x = gtu and the change of variables r = 
g/3{g). Recall that a{r) is the inverse of such mapping, furthermore, notice that dr = 
^{l + ^'{g))dg. This yields 

a(r) dr f , , f ip(z + a(r)u) , ( a(r)^ \ , 
^' ' ilj{roj)du ^ \\ dn,. 



Iq 1 + i!}'{a{r)) Jg2 J^± \z + a{r)uj\ \ \z + a{r)uj 

Turning back to cartesian coordinates a; = we obtain the desired expression 



^lj{x)TBi(p)ix) dx, 

with r B given by . □ 
The above result leads to a Carleman-like expression for ^: 

Corollary 4.2 (Carleman representation). Let e{-) be a restitution coefficient satisfying 
Assumptions ^J] and let 

B{u,a) = <^{\u\)b{u ■ a) 

be a collision kernel satisfying ( |2.6| ). Then, for any velocity distribution functions f, g one 
has 

QBeU.9){v)= [ fiz)[{t^oTBoQg]iv)dz 
Jr» 

where {t^ip] (x) = il){v — x) for any f , x G M'^ and test-function ip. 
Proof. The proof readily follows from the Lemma ^71] and the identity 

I QB,eU.9){v)i^{v)^v = \ I f{v)g{v-umu\)S_{Uij){u)dvdu (4.4) 

valid for any test-function ^. □ 
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4.2. Convolution-like estimates for ^. General convolution-like estimates are ob- 
tained in [|3|, Theorem 1] for non-constant restitution coefficient. Such estimates are given 
in LP with r] ^ and, for the applications we have in mind, we need to extend some 
of them to ?7 ^ 0. This can be done using the method developed in [ ]2T| ] (see also [|T4|]) 
together with the estimates of Q 

Theorem 4.3. Assume that the collision kernel B{u, a) = ^{\u\)b{u ■ a) satisfies ( |2.6D 



and ^{■) G L"^;^ for some /c G M. /n addition, assume that e(-) fulfills Assumption |2J 
Then, for any 1 ^ p ^ oo and ?7 G M, there exists Cn,p,k{B) > such that 

\\QUf.9)\\,.^C,MB) ll/IUUm. Il^ll-^... 
where the constant p ^(5) is given by: 

Cn,p,k{B) = Ck,^,p l{v,P,b) WHl^^ (4.5) 

— k 

with a constant Ck,r],p > depending only on k, r] and p. Furthermore, the dependence on 
the angular kernel is given by 

i{v,pM= / ^ Ks)ds> (4.6) 



where Ijp+l/p' = 1 andr]^ is the positive part of r]. Similarly, there exists C^^p^fc(-B) > 
such that ^ 

||Q5,e(/,^)IL.^c,,,,(5)||,IL.^^,^^^^, II/IU.^^ 

where the constant Cri,p,k{B) is given by 

Cn,p,k{B) =Ck,ri,pl{v,P,b) ||$||ioo (4.7) 

— k 

for some constant Ck,ri,p > depending only on k, rj and p. The dependence on the angular 
kernel is given by 

1 3+1+ 
7(r7, = + (1 - /3o)' ^ b{s) ds (4.8) 

where 1/p + 1/p' = 1 and /3o = /3(0) = 

Proof Fix 1 ^ p ^ oo and G M and use the convention 1/p' + 1/p = 1. By duality, 

||2^,e(/,^7)|L. = sup| / QlM.9)ivmv)dv : 
Using O, 

/" QBe(/>^?)WMdt^= / f{v)g{v~u)r4t,ij){u)dvdu 

with 

r-(V^)(M) = $(|m|)5_(^)(m), and t^V^(x) = V^(f - x), 



^Notice that the constants j{ri,p, b) and j{t],p, b) given by ( [4.6[ ) and (4.8) are not finite for arbitrary 
angular kernel b. It is implicitly assumed that the Theorem applies for the range of parameters leading to 



finite constants (see also Remark 
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with defined in equation With the notation of [Q], one recognizes that 5_ (h) 
V{h, 1), thus, applying [Q, Theorem 5] with q = oo and a = —i]. 



with j{ri,p, b) given by (p~6|). Notice that, with respect to IQ], we used the weight (t>)^ 
instead of this is the reason to introduce 77+ in our definition of 7(77,^,6). As a 
consequence. 



— Tj — k 



(4.9) 



Now, 



|/(t;)|dt;( / \g{u)\[{t,oT.otM{u)du 

'R3 



|/(t;)|||(t,or-ot. 



Using the inequality ||t>|Lp' ^ 2l^l/2(t;)NI||/i|| , for any s G R and v. 



dv. 



Jus 



^ 2\v+k\ + \v\/2^^^^p^ 5)||$|| 



— k~7] 



.1 dv 



dv 



ri + k 



which proves the first part of the Theorem. To prove the second part, observe that 



QlM,9)iv)^ljiv)dv= / fiv-u)g{v)r+Miu)dvdu, 

where T+{ip){u) = and S+ defined in (|4.ip. Using the notation of we 

identify S^{h) = V{1, h). Thus, applying Theorem 5] with p = 00 and a = —7], 



\h\ 



where ^{r],p, b) given by (p~8l). One concludes as above, interchanging the roles of / and 
9. □ 

Remark 4.4. The constants 'y{ri,p, b) and ■y{ri,p, b) are not finite for arbitrary be- 
cause of the possible singularity at s = ±1. However, if one assumes, as in [[191], that 
the angular kernel b{-) vanishes in the vicinity of s = 1 then '~f{r],p,b) < 00 for any 
1 ^ p ^ 00 and ?7 G M. This is an additional difficulty of the inelastic regime that is 
overcome in the elastic case using symmetry, i.e., defining b in half the domain. The care- 
ful reader will also notice that the constants given in the theorem are independent ofe{-) 
except for ^{rj,p, b) which depends only on the value e(0). Finally, we mention that, for 
hard-sphere interactions, i.e. b = one has 'y{r],p, b) < 00 <^==^ liv^P^ b) < 00 <^=^ 



1 ^p < 



3+r?+ 
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4.3. Sobolev regularity for smooth collision kernel. For this section we assume $( 
and b(-) smooth and compactly supported 



$ gc:o°°(m^\{o}), bec^{-i,i)- 

Denote by QB,e the associated collision operator defined by (I 



(4.10) 



Lemma 4.5. Assume that e(-) satisfies Assumptions |Z7| with e(-) G C'"(0, oo) for some 
integer m G N. Then, under assumption ( [4. 1(J| ) on the collision kernel, for any ^ s ^ m, 
there exists C = C{s, B, e) such that 



\TBif)\\Hs-,.^Cis,B,e) 



V/ G H' 



where Tb is the operator defined in Lemma 4J_. The constant C{s, B, e) depends only on 
s, on the collision kernel B and the restitution coefficient e(-). More precisely, C{s, B, e) 
depends on e(-) through the L°° norm of the derivatives D^e{-) (k = 1, . . . ,m) over some 
compact interval bounded away from zero depending only on B. 

We postpone the proof of Lemma ^ and first prove its important consequence. 



Theorem 4.6. Let B{u, a) = ■ a) be a collision kernel satisfying ( |4.10D and 

e(-) satisfying Assumption |2. 1\. In addition, assume that e{-) G ^"(0, oo) for some integer 
m G N. Then, for any ^ s ^ m, 



\QteU,9)\\^.^.^C{s,B,e) \\g\\HA\f\W 
with constant C{s, B, e) given in Lemma [t.5[ 



Proof. Let T [Q^ g)] (^) denote the Fourier transform of ^{f, g). According to 
Corollary I 



^[QBAf^9)]{0= / f{v)J'[{t.orBoQg]{Odv. 
To simplify notation set G{v,^) = T \{t„ o Tb o i^j) S'] (0- Thus, 



(0 



2(s+l) 



f{y)G{v,i) dt; 



(4.11) 



\f{v)\\G{v,i)\\if^^^^^&i&v. 
Since G{v, i) = 7 \{t, o o t,) g\ (0, 

/ \G{y, di = ° o Q g\\%^, ^ G{s, B, ef Wgf^. 

JR3 



For this inequality we used Lemma ^3] and the fact that the translation operator has 
norm one in any Sobolev space. Hence, estimate ( |4.1 1| ) yields the desired estimate. □ 
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Proof of Lemma The proof of the regularity property of Tb can be obtained following 
the lines of the one for the elastic Boltzmann operator [|2l|]. Indeed, note that 

r^(/)(r,a;) : = V BU){o^'\r),uj) = TB{f){rP{r),co) 
B{z + ru, r)ip{ruj + z) dvr^. 



Assumption ( [4.1Q[ ) implies that there exists 5 > such that b{x) = for |a; ± 1| ^5 and 
{|z| ; z E Supp($)} C (a, M) for some positive constants < a < M. Then, by virtue 
of iQ,B{z + 



ruj.r] 



for any r > 0, w G § and z E uj provided that \z\ > ■ 

For \z^ ^ o'^^ has |2; + rwp ^ 2r'^ /5, thus, + rci;,r) = if r < ^j5a^/2. 

Putting these together we conclude that 



B{z + rw, r) = Vr ^ / := (^'^6a^/2, M j , w G and any ^±tj. 
In particular, V B{f){r,uj) = for any r ^ / independently of /. Define 



(4.12) 



Bo{z,g) ■.= l±fi^/3\g)B{z,g) 



H\z\)b{l-2 



e 



and denote ro(/) the associated operator. 



ro(/)(r,u;) := 



Bq{z + rcj, r)(y9(rcj + z^ dvr^. 



Then, i3o does not depend on the restitution coefficient e(-) and Fq is exactly of the form 
of the operator T studied in [|21], Theorem 3.1]. Therefore, arguing as in op. cit., for any 
s ^ 0, there is an explicit constant Co = Co(s, $, 6) such that 



ro(/) 



^Co(s,<l>,6) 



V/ G H'. 



Setting 



^ 0, 



(4.13) 



(4.14) 



one observes that Ge is a C"^ function over / whose derivatives D'^Ge are bounded over / 
for any k ^ m and 

rB{f){r,co) = G,{r)xi{r)I\,{f){r,co). 



Here xi is the characteristic function of / = \^^/6a?J2, Mj (see Eq. ( [4-.12D ). Therefore, 
for any ^ s ^ m, there exists some constant G = Go{s, b, e) such that 



Hs+i 



^Gois,B,e) 



V/ G H' 



(4.15) 



where the constant Go{s, B, e) can be chosen as 



Co(s, 5, e) = Co(s, $, 6) max \\D''Ge\\L'^ii). 

k=0,...,s 



(4.16) 
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From estimate ( |4.15D we deduce Lemma with the following argument. Assume first 
s = /c ^ 1 is an integer. Using polar coordinates 

/»oo I* 

\\^Bif)\\l. = Y. F,ig)g'dg \diTBif)ig,u)\' du 

where, for any |j| ^ k, the function Fj(g) can be written as 

F,{g) = P,(^«(f?), • • • + ^'^'\g))-''^ . (4.17) 

Here Pj{yi, . . . ,yj) is a suitable polynomial, uj G N and 'd^^^ denotes the p-th deriva- 
tive of ■{?(■). Since ^9 G ^"(0, oo) and / is a compact interval away from zero, one has 
sup^gj Fj{g) = Ck < oo for any \j\ ^ k. Thus 

\\TBif)\\H^^Ck\\TBif)\\H^ (4.18) 

where Ck is an explicit constant involving the L°° norm of the first k-th order derivatives 
of «(■) on /. This proves that the conclusion of the Lemma p3| holds true for any integer 
s ^ m and we deduce the general case using interpolation. □ 

Remark 4.7. It is important, for our subsequent analysis, to obtain a precise expression 
for the constant C{s,B,e). For instance, in the case in which e(-) G C^(0, oo), one 
obtains that 

C{l,B,e) ^Co{l,B,e) supFi(f?) 

eel 

where Fi is of the form ( |4.17[ ) with I defined in ( [4.12| ). Note that Co{l,B,e) and Ge{g) are 
given by ( |4.16| ) and ( |4.14D respectively. In particular, under Assumption Ge{g) ^ 4p 
for large g and Ge{g) — g/2for ^ ~ 0. 

Arguing as in [ pT] , Corollary 3.2] we translate the gain of regularity obtained in Theo- 
rem in gain of integrability. 

Corollary 4.8. Let B{u, a) = ^{\u\)b{u ■ a) be a collision kernel satisfying ( |4.10D and 
e(-) G C^(0, oo) satisfying Assumption |2. I\ Then, for any 1 < p < oo 

\\QteU\9)\\,r> ^ C{p,B,e) {\\gU. \\fh^ + \\gh^ ||/|U0 

where the constant C{p,B,e) depends on B and e through the constant C{l,B,e) of 
Theorem 4^. The exponent q < p is given by 

JTTp (4.19) 
p/3 if pE [6, oo). 

4.4. Regularity and integrability for hard-spheres. We consider in this section the 
case of hard-spheres collision kernel 



B(u,a) 



\u\ 
Att' 

Such a collision kernel does not enjoy the regularity properties assumed in the previous 
section. This does not present a problem since the dependence of the constant on the 
collision kernel B permits to adapt the method developed in [ pi] ] for the elastic case. We 
need some supplementary assumptions on the restitution coefficient e(-). 
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Assumptions 4.9. In addition to the Assumptions |Z7|, suppose that e{ ) G C^(0, oo) and 
that there exists /c G M such that 

e'(r) = 0{r^) when r — > oo, 

where e'(-) denotes the derivative ofe{-). 

The above assumptions imply that '&'{g) = 0{q^^^) for large q and ^'{q) ~ 1 when 
^ ^ 0. Recall that ^'{■) is the derivative of d{r) = re(r). 

Theorem 4.10. Assume that e(-) satisfies Assumptions \4.'^ For any p G [1, 3) there 
exist K > 0, G (0, 1) and a constant Ce > depending only on p and the restitution 
coefficient e(-) such that, for any 6 > 



[ Qtif, f) dv ^ CJ-^ 11/11 11/11?;-'^ + 5 



i/p 



Remark 4.11. The restriction p E [1,3) is the major difference with respect to the classi- 
cal case [ pi] , Theorem 3.1]. The reason is that in the inelastic regime the lack of symmetry 
does not permit to switch the roles ofv' and v'^, therefore, general b has to be defined in 
the full interval [—1, 1]. 

Proof. We follow the same lines presented in [|T]] and subsequently used in [[T^]. We 
present the argument for convenience. Fix p E [1,3) and let 6 : R — t- IR+ be an even C°° 
function with compact support in (—1, 1) and J^^ 6(s) ds = 1. In the same way, consider 
a radial C°° function E : — )■ M with support in the ball B{0, 1) and Jj^g S(f ) dv = 1. 



Define the mollifications S„ 

'^Sn = S„ * (I • \xaJ and 6c 



n^E{nv) and 6m(s) 



4lT 



L m. ' rn, J 



collision kernel. Here we have defined the set 



V G 





'2 






-,n 


j- n 




n 




a) = $ 


s.{\u 


)bsm{u- 



:= m&{ms) for m, n ^ 1. Thus, 
are smooth mollifications of the 



Consider the smooth collision kernel 

and observe that 

supp (<l>5„) C |i ^ ^ n + 1 [> and supp (bsj C 
Define naturally 



n 



-1 + 



1 - 



m 



m 



Here<l>/j„(|?i|) 



BsR^,AH,u- a) 

Brr^,„{\u\,u- a) 
|?i|-$5„(|u|)and 6ij„ 



"^RniluDbs^Xu-a) and 
^RM)bR^iu-a). 



(u-a 



4lT 



— bs^ {u ■ a) are the remainder parts. 



Thus, one splits in four parts using obvious notation. 



Q 



Brr„ 



HAFF'S LAW FOR VISCOESLASTIC HARD SPHERES 



29 



Since Bs^„{\u\,u-a) fulfills ( |4.1(]| ) one deduces from Corollary p~8|that there is a constant 
C(m, n) such that 



Q 



+ 



LP 



Li 



for g < p given by ( [4.19| ). A simple application of Holder's inequality yields 

2^s_,e(/>/) < C{m,n) ll/IU, ll/IU. ||/||^;\ (4.20) 

Recall from Corollary ^]8] that C(m, n) depends on m and n through the constant C(l, Bs^„, < 



in Theorem pTq . Moreover, according to Remark pk7| , one sees that 

C(l, e) ^ Co(l, <^s„,bsj max ||D'=Ge||L-(/) sup Fi{q) 

fc=0,l 

where Co(s, $, 6) is the constant appearing in ( |4.13p , Ge{-) is given by ( |4.14| ), and Fi is 
of the form ( |4J7| ). The interval / = is defined in with 5 = l/m, M = n + 1 

and a = 1/n 

,n + 1 



2m'n? 



That Co(l, $5„, 65^) depends on m and n in a polynomial way follows as in [PT]]. More- 
over, from the properties of Ge given in Remark and the fact that Fi{q) is a rational 
function in one deduces from Assumption ^]9| and the above expression of / that 

there exist a, 6 > such that 



C{m,n) = 0{m^ n^) as m, n — )■ 00. 
Now, applying Theorem [4 . 3| with k = 1 and rj = — we get 



(4.21) 



Q 



.if J) 



(/,/) „^^^o("^,^ 



1/p 



where eo{m,n) = C_i/p/,p,i(55/j„ „) + C_i/p/,p,i(5Ri?„ „) forany m,n ^ 1. In particu- 
lar, using the expression of the above constants in (p3|), there exists a constant c > such 
that eo{m,n) ^ c'y{—l/p',p,bji^) =: e(m) for any m, n ^ 1. Then there exists some 
r > such that 

e{m) = 0{m~'^) as m — )■ cxd. (4.22) 

Indeed, since 1 ^ p < 3, one sees from (^3j) that 7(— bji^) ^ C'Ufei?^ ||l9(§2) for 
any q such that 1 < g' < 2p' /3. Thus, one can choose a regularizing function 6 so that the 
L^(S^)-norm of bn^ decays algebraically to zero as m grows. Using the above estimate 
with 1] = we get 



It remains only to estimate 



m 



LI 



(4.23) 



I : = 



Q 



R3 



{f,f)r~'dv. 



One notes that 



<l>RA\v-vA)^Cn-mv\' + \v.\') 
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for some C > 0. Thus, 



I ^ Cn-^ 



Define 



Ii : 



f{v)f{vM'dvdv^ r-\v')bsju-a)dcr, and 
fiv)fiv^)\v,\'dvdv^ [ r-\v')bsju-a)da. 



Observe that Ii can be written as 

Ii = 

where 



QtAFJ){v)ij{v)dv 



F{v) = \v\^f{v), tP{v) = r-\v) e LP'{R^) 

with the collision kernel Bra{\u\,u - a) = hs^iu ■ a). Applying Theorem ^ with rj 
k = gives 

^ CoMBn.) WFhAlfhrmiL.' < Co,p,o(5™)||/||i 
where Co,p,o(-Bm) is defined by (|4.5D. Now, with the same notation. 



LP 



QtAf^F){v)ij{v)dv, 



therefore, applying Theorem with r] = and k = —2 yields 

h ^ Co,pMBm) II/IIlI lli^lk^ Mlp' < Co,p,-2(5„ 

Combining the two estimates for Ii and I2, 

C{m) 



p 

LP- 



I ^ 



n 



L\ 11/ II LP 



where C{m) = Co,p,o(-Bm) + Co,p,-2(-Bm)- The support of hs^{s) lies to a positive 
distance, of order 1/m, from s = 1. Then, we use the expression ( |4.5| ) to conclude that 



C{m) ^ m 27 as m — i- 00. 
Estimates ( F^ , d^IOl ) and gives 



(4.24) 



C{m) 



LP- 



Using the polynomial bounds ( [4-.2ip , ( |4.22| ) and ( |4.24| ) this leads to the result as in [fT^]. 

□ 
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Remark 4.12. Assumption \l. 9| allows to present the explicit dependence of the constants 
with respect to 5 > 0. This dependence will be crucial in the proof of Haff's law in Section 
5. Note that the constant Ce in Theorem }^.10\ depends on the regularity of the restitution 
coefficient away from zero. 

Corollary 4.13. Assume that e(-) satisfies Assumption p79[ For any p E [1, 3) there exist 
K > 0, G (0, 1) and a constant Ce > depending only on p and the restitution coefficient 
e(-) such that, for any 6 > 

[ Qtig,9)g'-\vrdv^CJ'^\\g\\ir''\\9C-'h5\^^^^ \\g\\l. , 0. 
The constant Ce is provided by Theorem \1.1C\ 

Proof Fix 5( ^ 0, ^ and set f{v) = g{v){v)'^. Note that {v')'^ ^ (v)"^ (v^)'' for any 
V, Vi, G M'^, then, using the weak formulation of 

/ Qt{9.9)f-\vy'^<lv= ! {vrQt{9,9)F^'dv^ [ Qt{fJ)f''-'dv. 

Jr^ JR3 J-jg3 

Conclude with Theorem |4.10| . □ 

The following result applies to the rescaled solutions g{T, w). Its importance lies in that 
the estimate is uniform in the rescaled time r. 



Corollary 4.14. Assume that e(-) satisfies Assumption For any r ^ 0, let be 

the restitution coefficient defined by ( |2.13D and let Qe^ (/, /) be the associated collision 



operator Assume that V{C{t)) is continuous and goes to infinity as t ^ oo. For any 
p G [1, 3) there exist k > 0, 9 E {0, 1) and K > all independent of t such that, for any 

[ Qt{9,g)g''\wrdw^K6--\\g\\lt^'\\gr^^^^^ \\g\\l, , Vr/ ^ 0. 



Proof. From Corollary [4.13| , for any r ^ there exists K{t) = Cz^ for which the above 



inequality holds. It suffices to prove that K = snp^^Q K{t) < oo. Recall that K{t) 
depends on r through the restitution coefficient Cr, more precisely, Ce^ depends on the 
L°° norm of the derivatives D'^e^(-), k = 0, 1, over some compact interval of (0, oo) 
bounded away from zero (independent of r). Now, for any r ^ 0, 



with p{t) = V{C{t)). Since p.^^{r) is continuous and goes to zero as r goes to oo, one 
concludes that all the norms of D''e^(-) remain uniformly bounded with respect to r. 
The same holds for K{t). □ 

5. Generalized Haff's law continued 
5.1. Proof of Haff's law. In this section we prove the second part of Haff's law es- 



tablishing the lower bound of the temperature (p.l2D. Recall that, from Theorem P]7| it 



suffices to prove ( p.lO[ ). As explained in Section ^ this is done using suitable U' estimates 
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in the self-similar variables. In this section, the restitution coefficient fulfills Assump- 
tions ^J] and ^]9| and the collision kernel is that of hard-spheres interactions. Recall that 
the rescaled function g{T, w) is solution to the Boltzmann equation in rescaled variables 
CT ) 

drg{T,w) + ^{T)V^-{wg{T,w)) = Q^^{g,g){T,w) r > 0. (5.1) 
The restitution coefficient and the time-depending mapping ^(r) are given by ( |3.19| ). 

Proposition 5.1. Assume that e(-) fulfills Assumptions with 7 > and pT^ Let fo 
satisfying ( pl8| ) with /o G fl LP{M.^)for some 1 < p < 3. Let g{T,-) be the solution to 
the rescaled equation ( |5.1| ) with initial datum g{0, w) = fo{w). Then, there exist Co > 
and /to > such that 

Mr)\\L.^Co{l + TY' Vr^O. (5.2) 
Consequently, there exist Ci > and ki > such that 



g{T,w)\w\^ dw ^ Ci(l + t)~ 



Vr > 0. 



(5.3) 



Proof. The proof relies on Corollary [4.14 Multiply ( |5.1D by (7^ ^ and integrate over 
to obtain 



ld||^(r)||i,^3/ ^ 



dr 



p 



^r) \\g\ 



T 



[ Qti9,9)9'''dw- [ Q-ig,g)g^-'dw. (5.4) 



Vr > 0. 



(5.5) 



From Jensen's equality, one has 

/ Q~ig,g)g''-'dw^ [ gPiT,w)\w\dw 

According to Corollary |4 . 1 4| there exist k > 0, 9 E (0, 1) and a constant K > that does 
not depend on r such that 



> 0. 



JR3 1/p 

From conservation of mass ||5'(t)|| = 1, furthermore, M2 := sup^^Q \\g{r) ||^i < 00 from 
(|3Tl7|) . Thus, using (jsj) and (F^b , 



d|k(r) 



I LP 



^ pirr«||^(r)ri^-^) + pM2 51|^7(r)r^. - Mr)lk(r)r^p (5.6) 

QT i/p i/p 

where /i(r) = min (1, 3{p — l)^(r)). Since ,^(r) — t- as r — t- 00 for 7 > 0, there exists 
To > such that 

/i(r) = 3{p - l)^(r) = _y 7 for any ^ ^ ^0- 



7r + 1 + 7 



Choosing 5 = ^{t)/ {PM2) in (K^ we get 



d|k(r) 
dr 



I LP 



^pir(pM2)>(r)-«||^(r)||^(;-^) ^ C'(7r+l+7)1|^7(r)||ii^-'') Vr ^ 



To 
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for some positive constant C > 0. Integrating the above estimate, we conclude the exis- 
tence of some constant Cq > such that 

\\9{r)\\%^Co{lT + l + ^)'^ Vr^ro, 
and ( |5.2| ) readily follows. 

Regarding estimate (^31) note that for any R > 0, 



@(r) 



g{T,w)\w\'^ dw + 



|«>|>-R 
2 



^\w\>R 

From Holder's inequality, 



gij, w)\w\ dw 



Vr ^ 0, 



g{T, 'w)\w\ dw ^ 



/4 y/p' 1 1 

( -vri?^ I IkfTlllLp with the convention — I = 1. 

V3 / V V' 



Therefore, using (p^, there exists a positive constant C > independent of R such that 

0(r) ^ R^ (\ -CR^I^'il + r)"«) Vi? > 0, Vr ^ 0. 
Pick i? = R{t) > such that Ci?3/p'(i + t-)«o = 1/2, then 



@(r) > ^-R^{t\ 



1 



2 V2C(1 + r) 



p'/3 



Vr > 0, 



which gives ( |5l3| ) with k\ = p'kq/?,. 



□ 



The generalized Haff 's law is a consequence Theorem ^]7] and Proposition |5TT 



Theorem 5.2. Le? /o ^ satisfy the conditions given by ( ^l8| ) wzY/z /o G L*'°(M^) /or 
5ome 1 < po < 00. /n addition, assume that e{ ) fulfills Assumptions 3J_ and |4. 9[ T/zen, 
f/ze solution f(t, v) to the associated Boltzmann equation ( |2.7[ ) satisfies the generalized 
Haff's law 

c(l + t)~^ ^^(t) + t^O (5.7) 

where c, C are positive constants depending only on e(-) and S{0). 

Proof. The upper bound in ( |5.7D has already been obtained in Theorem The proof 



of the lower bound is a straightforward consequence of Theorem [377| and Proposition |5.1|. 
Indeed, notice that if /o G L^(R^) fl Lp°(M.^) for some 1 < < 00, using interpolation, 
we may assume without loss of generally that po G (1, 3). Recall that for 7 > 0, 

£{t) = V~'{t)&{r{t)) 

where V{t) = (1 + t) 1+^ and r(t) is given by ( P.18D . Since ©(•) decays at least alge- 
braically (p3|), one recognizes that there exists some constant a > such that £{t) ^ 
a (1 + with /i = ^±2^ with ki being the rate in (|5.3|). The result follows from 

^ □ 



Theorem |3/7|. The proof for 7 = is identical. 
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Example 5.3. For constant restitution coefficient 7 = 0, we recover the classical Haff's 
law o/[|T5l] proved recently in [fT^]; 

Example 5.4. For viscoelastic hard-spheres given in Example p.^ one has 7 = 1/5. 
Thus, Theorem ^^provides the first rigorous justification of the cooling rate conjectured 
in 

c(i + t)-^/=^ ^C(l + ^)-^/^ t^o. 

Remark 5.5. Theorem shows that the decay of the temperature is governed by the 
behavior of the restitution coefficient e(r) for small impact. The cooling of the gases is 
slower for larger 7. 

From the explicit rate of cooling of the temperature, one deduces the algebraic decay 
of any moments of the solution to (|2l7|). Under the assumptions of the above Theorem |5^ 
the p— moment nipit) defined in ( p.4p satisfies 

Cp{l + t)-^ ^£{tf i^mp{t) ^Cp£{tY ^Cp{l + ty^ , t^O. (5.8) 

The positive constants Cp, Cp, Cp depend on p, mp{0), £(0) and e(-). The lower bound 
is a direct consequence of Jensen's inequality and (^Jh while the upper bound has been 



established in Theorem 3.7 



5.2. Application: Propagation of Lebesgue norms. We complement Proposition |57T| 
by proving the propagation of L^-norms in the range 1 ^ p < 3 for the solution g(T, w) 
satisfying the rescaled equation (|5.1[). Thus, the method introduced in the elastic case 
pT} ] and later used in [ p^ for constant restitution coefficient is extended to the case of a 



variable restitution coefficient satisfying Assumptions ^7T] and 



Lemma 5.6. Assume that the initial /o ^ satisfies the conditions given by (|2.8D with 
/o € L^(M^) for some 1 < p < 00 and let g{T, ■) be the solution to the rescaled equation 
( |5.1D with initial datum (7(0, w) = fo{w). Then, there exists a constant z/q > such that 

/ g{T,w^)\w ~ w^\dw^ ^ ma.x{i'Q,\w\} ^ —{w), Vty G M'^, r > 0. 
Jr3 2 

In particular, 

[ 9'-'Q:{9,9)dw^^ [ g^{T,w)il + \w\'Y/'dw='^MT)\\l, . 

Jm.3 I J^3 I 1/P 

Proof. The proof is a simple consequence of 

©mill := inf / (7(r, dw > 0. 
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Indeed, since /q G the propagation of p-moments in the rescaled variables implies 
supt^o < Then, for i? > large enough 

/ g{T,w)\w\'^ dw = / g{T,w)\w\'^dw— / g{r,w)\w\'^ dw 

J{\w\^R} Jr^ J{\w\^R} 

> Qmin - 4 sup ||^(r)||^i ^ 0min/2 > 0. 

We conclude that, 

5((r,t(;)|w| dw ^ 4 [ g{T,w)\w\'^dw^^:^=:uo>0. 

Using this observation and Jensen's inequality we obtain the result. □ 

Theorem 5.7. Assume the variable restitution coefficient e{-) satisfy Assumptions [3. 1\ and 
^^for some positive 7 > 0. Assume that /o ^ satisfies with /o G Ll^-^_^^-^r\LP(R^) 
for some 1 ^ p < 3 and rj ^ 0. Then, the rescaled solution g{T,-) to ( ]5.1D with initial 
datum g{0, w) = fo{w) satisfies 



sup \\gir)\\^p < 00. 

In particular, 



sup \v{t)-'/p' > = sup \\g{T)\\^, < 00. 

Recall that V{t) = (l + t)^. 

Proof. Multiplying equation [STTl by gP~^{T, w) {w)^'^ and integrating over M'^ yields 

idlkMlli. 



- + ^(^-l)ar)\\9\\l.= [ Qt{9,9)9 



{wY'' dw- 



p dr 



Using Lemma one has 

/ Q-{9,9)9'-'{wr dw^^'-^U 
Jr3 ^ 



ri + l/p 



Moreover, = sup^^o < 00 by virtue of the propagation of moments in 

4C 



self-similar variables (|5^). Applying Corollary [4.14] with 5 



-^lk(r)||^ + ^IKr)r,. 
p dr 4 ^v+i/p 



^^IKr)r,r^)+e(r)(ry-^)||^(r)r,. Vr > (5.9) 

for some uniform constant K. Since 7 > 0, the mapping ^{r) decreases toward zero, thus 
(^3* l^^d^ to the result. □ 
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Remark 5.8. We refer to [19, Theorem 1.3] /or a proof of the case 7 = 0. Furthermore, 
additional pointwise estimates allow to extend the above result to p ^ 3 assuming higher 
moments for /q. We refer the reader to ^ for further similar estimates. 

6. High-energy tails for the self-similar solution 

We finalize this work studying the high-energy tails of f{t, v) of the solution to (|1.3|). 
For models with variable restitution coefficient the high energy tail is dynamic since gas 
changes its behavior during the cooling process. This is noted with a dynamic rate in the 
tail. Here again, we shall deal with the generalized hard-spheres collision kernel 

B{u, cr) = \u\b{u ■ a) 



where b{-) satisfies ([2.6D. We argue in the self-similar variables, thus it is convenient to 
define the rescaled p-moments 

™-pi'^) = / g{T,w) \w\'^^ dw, p^O. 

Notice that ( |5^ readily translates into 

Cp ^ mp(r) ^ Cp for r ^ 0. (6.1) 

The following Theorem generalizes [[T9|, Proposition 3.1] to the case of a variable restitu- 
tion coefficient. 

Theorem 6.1 (L^ -exponential tails Theorem). Let B{u,a) = \u\b(u ■ cr) satisfy ( |2^ ) 
with b G L'^ iE>^) for some q > 1. Assume that e{-) and fp fulfill Assumptions [j. i| and ( |2.g| ) 
respectively. Furthermore, assume that there exists rg > such that 

fo{v) exp {ro\v\) dv < oo. 

Let g{T, w) be the rescaled solution defined by ( |2.1Q| ). Then, there exists some r ^ tq such 
that 

sup / g{T,w)exp{r\w\) dw < oo. (6.2) 

T^O Jr3 

Consequently, 

sup / /(t, f ) exp (rV"(t)|f I) dw < oo. (6.3) 

Proof. The method of proof is carefully documented in |^] . We sketch the proof divid- 
ing the argument in 5 steps. 

Step L Note that formally 

« oo J. 

/ g{T, w) exp (r|w|') dw = 22 Tj-msA:/2(r), 

"^^^ A:=0 



for any r > and any s > 0. Hence, the summability of the integral is described by the 
behavior ( 
moments 

^.Jt^ ■ = 

T{ap + b) 



behavior of the functions — ^^^p-^- This motivates the introduction of the renormalized 

^p(^) •= T^fl\ with a = 2/s, 
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where r(-) denotes the Gamma function. We shall prove that the series converges for 
some r < To and with s = 1 (i.e. a = 2). To do so, it is enough to prove that, for some 
6 < 1 and Q > large enough, one has Zp{T) ^ for any p ^ 1 and any r ^ 0. 



Step 2. Recall that, according to Lemma the estimates of Proposition |2J| are inde- 
pendent of the restitution coefficient e(-). In particular, they hold for the time-dependent 
collision operator Q^^ providing bounds which are uniform with respect to r. Specifically, 

/ QeA9,9){T,w)\w\^Pdw ^ -{1 ~ Kp)mp+i/2{T) + K,pSp{T), Vt ^ 

where Kjp IS the constant introduced in Lemma and 



Mr) 



( I ) (™fc+i/2(^) ^P-k{r) + mfc(r) mp_k+i/2{r)) . 



Step 3. An important simplification, first observed in [§], consists in noticing that the 



term Sp satisfies 

5p(r) ^ A T{ap + a/2 + 2h) Zp^r) for a ^ 1, b > 0, 
where A = A{a,b) > does not depend on p and 

Zp{r) = max {zfc+i/2(r) Zp^kir), Zkir) Zp^k+i/2ir)} . 

With such an estimate, the rather involved term Sp is more tractable. 



Step 4. Using the above steps and the evolution problem ( |5.1D satisfied by the rescaled 
solution g, we check that 

dm 

-^(r) + (1 - Kp)nip+y2{r) ^ Kp r{ap + a/2 + 2b)Zp{T) + 2pe(r)mp(r) 
where we used the fact that 

/ \w\'^''Vw ■ {wg{T,w)) dw = —2pmp{T). 

Using the asymptotic formula 

l.m = 1. 

p-s>oo L [p + S) 

the fact that ^(r) ^ 1 and Kp ~ l/p^^''' for large p, one concludes that there are constants 
Cj > (i = 1, 2) and Po > I sufficiently large so that 

^(r) + ci p'^/2^;+V2p(^) ^ C2 p-Z^+'-yi' Zp{r) + 2p Zp{r) Wr ^ 0, p ^ p,. 

We also used that mp_|_i/2(T) ^ mp^^^^^(r) for any r ^ thanks to Jensen's inequality. 

Final step. We claim that if we choose a = 2 and < b < it is possible to find 
Q > large enough so that mp(r) ^ Q^. Indeed, let po and Q < oo such that 

—Po'^^'^ ^ t:, and I max sup ZA;(r), Qo, 1 
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where Qq is a constant such that Zp{0) ^ Q^. This constant exists by the exponential 
integrability assumption on the initial datum. Moreover, since moments of g are uni- 
formly propagated, the existence of such finite Q is guaranteed. Arguing by induction and 
standard comparison of ODE's, one proves that Upir) := Q'^ satisfies for p ^ po 

therefore, yp{r) ^ Zp{T) for any p ^ po. Since this is trivially true for p < powe obtain 
that 

mp{T) ^T{2p + b)QP, Vp^l,r^O. 
From Step 1, this is enough to prove the Theorem. □ 

Example 6.2. For viscoelastic hard-spheres V {t) = (1 + t)^/^. Therefore, 

/ fo{v) exp (ro|f I) df < oo =^ sup / f{t, v) exp (r(l + t)^''^|t>|) dv < oo 

for some r < vq. In particular, using the terminology o/ [Q], fit,v) has a (dynamic) 
exponential tail of order 1. 

Appendix A: Viscoelastic hard-spheres 

In this Appendix we prove that Assumptions ^71] are met by the restitution coefficient 
e(-) associated to the so-called viscoelastic hard-spheres as derived in p2| ] (see also [ |T^ , 
Chapter 4]). In fact, we prove a more general result for the hard-spheres collision kernel 

B(u,a) = — WeR^ae^"^. 

4:71 



Recall that *e was defined in (|3.1| ) as 

= /""^ (1 - dz, a; > 0. 

Lemma A. 1. Assume that e(-) satisfies Assumption ^.l\ and that the mapping r ^ 
e(r) is decreasing. Then, the associated function defined in ( p. ID is strictly increasing 
and convex. 

Proof. Since e is decreasing, e'(r) ^ for any r ^ 0. Here e'(-) denotes the derivative of 
e(-). Define 

:= - / (l-e'^(z)) z^dz, x > 0. 
X Jo 

Note that ^'e(') is convex if and only if x^xx{x) — ^xi^) ^ for any x > where $3, 
and denote the first and second derivatives of $ respectively. A simple calculation 
shows that 

3 r 

x^xxix) - ^x{x) = -2x^e'{x)e{x) + — (1 - e'^{z))z^ dz, Vx > 0. 

^ Jo 

Since e'(x) ^ and e(-) G (0, 1] one concludes that x^xxix) — $^(0;) ^ for any x > 0. 

Similarly, since e'(-) ^ the mapping z ^ (1 — e^{z))z^ is nondecreasing, thus, 
^xix) > for any x > 0. This implies that *e(") is strictly increasing over (0, +00). □ 
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For the visco-elastic hard-spheres, as derived in p2[], the restitution coefficient e is 
solution of the equation 

e(r) + ar^/^e(r)3/^ = 1 Vr ^ (A.l) 

where a > is a constant depending on the material viscosity. It was proved in p. 
1006] that, on the basis of ( |A.1| ), Assumptions ^]l] are met. From equation ( JA.ID , one 
deduces that 

lim e(r) = 1, and e(r) ^ 1 — ar^^^ for r ~ 

which means that Assumption ^71] (i j is met. Furthermore, equation ( |A.1| ) also implies 
that e is continuously decreasing. According to Lemma A.|T], e(-) satisfy Assumptions |3.1|. 
Moreover, it is easy to deduce from ( |A.1D that Assumption is satisfied. 

Example A. 1. For monotone decreasing restitution coefficient introduced in Example 
|2Ji Assumptions^^^are also met by virtue of the above Lemma. In such a case, according 
to ( |2.2| ), the cooling of the temperature £{t) is 



Sit) = 0[{l + t) ^) as t^oo. 
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